Guaranteed and Prospective Galactic TeV Neutrino Sources 
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Recent observations, particularly from the HESS Collaboration, have revealed rich Galactic populations of 
TeV gamma-ray sources, including a collection unseen in other wavelengths. Many of these gamma-ray spectra 
are well measured up to ~ 10 TeV, where low statistics make observations by air Cerenkov telescopes difficult. 
To understand these mysterious sources, especially at much higher energies-where a cutoff should eventually 
appear-new techniques are needed. We point out the following: (1) For a number of sources, it is very likely 
that pions, and hence TeV neutrinos, are produced; (2) As a general point, neutrinos should be a better probe of 
the highest energies than gamma rays, due to increasing detector efficiency; and (3) For several specific sources, 
the detection prospects for km 3 neutrino telescopes are very good, ~ 1-10 events/year, with low atmospheric 
neutrino background rates above reasonable energy thresholds. Such signal rates, as small as they may seem, 
will allow neutrino telescopes to powerfully discriminate between models for the Galactic TeV sources, with 
important consequences for our understanding of cosmic-ray production. 

PACS numbers: 95.85.Ry, 95.85.Pw, 98.70.Rz, 98.35.-a 



I. INTRODUCTION 

The field of TeV gamma-ray astronomy is exploring en- 
ergy regimes that have been, until recently, out of reach to 
astrophysicists. Yet, even as the catalog of TeV sources 
continues to grow, it is still debated whether the observed 
gamma rays are produced leptonically, through the inverse 
Compton scattering of energetic electrons on ambient pho- 
tons (e~y — > ye~), or hadronically, though neutral pion de- 
cay (n° — > 77). Air Cerenkov telescopes, such as HESS 1 1], 
can measure a source spectrum with high precision for Ey ~ 
1 — 10 TeV. One might hope that as we probe higher en- 
ergy gamma rays, indications of their true origin would be 
revealed. However, at energies > 10 TeV, the difficulties of 
gamma-ray astronomy become more pronounced, due to the 
low statistics of the quickly declining signal spectra. 

It is well established that a distinctive feature of a pionic 
(hadronically-produced) gamma-ray spectrum is an accompa- 
nying flux of neutrinos |2j yl U, |5[ . These neutrinos originate 
from the decay of charged pions (7T+, 7Z~), which are pro- 
duced in approximately equal numbers with neutral pions in 
proton-proton scattering. Relative to gamma-ray telescopes, 
the new km 3 neutrino telescopes 0, 0] will have several ad- 
vantages that result in improved performance at the high- 
est energies. The rapidly falling atmospheric neutrino back- 
ground, rising neutrino-nucleon cross section (<J v n ~ E v ), an d 
increasing muon range, which effectively expands the (al- 
ready large) detector volume (R^ ~ InE^), all help to amplify 
the diminished flux at these energies. In fact, as the back- 
ground quickly becomes negligible in the TeV range, the de- 
tection of any high energy neutrinos from a source could sig- 
nificantly indicate a hadronic production mechanism. 

Neutrino telescopes have capabilities far beyond break- 



ing the degeneracy between leptonic and hadronic production 
models. Spectral features in the highest energy regime, espe- 
cially an expected cutoff (related to the maximum accelerated 
proton energies), shall not remain inaccessible to observation. 
Measurement of the energies of neutrino-induced muons and 
showers, which are related to the original charged pion en- 
ergy, can probe the source proton spectrum in a complemen- 
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FIG. 1 : Ranges over which gamma-ray and neutrino observations 
can measure a source spectrum (Ey ~ E n /2, E v ~ E n /A). Squares 
are HESS measurements of Vela Jr., fit by a power law (dashed line). 
Neutrino telescopes can probe higher energies to distinguish between 
pionic spectra with cutoffs (dot-dashed lines) by measuring the - 
induced muon and shower spectra. Muon rates of N(Eu > 1 TeV) 
~ 3 — 6 yr~' correspond to the dot-dashed lines (see Section HW 
illustrating the sensitivity of neutrino observations. 
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tary manner to gamma-ray observations, which effectively 
measure the neutral pion spectrum The sensitivities 

of these two independent approaches, including the regime 
where they coincide, are illustrated in Fig.Q] The ability to ac- 
curately measure neutrino-induced muon spectra greatly im- 
proves the prospects for detecting point sources, as the harder 
source spectra dominate the atmospheric background above 
~ 1 TeV. While muon tracks have better angular resolution 
(< 1°), neutrino showers (~ 10° in water) more faithfully 
trace the spectrum. Shower observations, which measure the 
V e and v t fluxes, when combined with muon data, also allow 
for the study of the ratio of neutrino flavors arriving from a 
source 

Considering the latest observations of Galactic TeV 
sources, we calculate the corresponding spectra of detectable 
neutrino-induced muons and showers, assuming only that the 
observed gamma-ray spectra are pionic, for a range of pos- 
sible high energy cutoffs in the spectra. Relative to analyses 
in which the total numbers of signal and background events 
are counted (e.g., Ref. 112]), a maximum likelihood analysis 
would have much more power. For example, below ~ 1 TeV, 
a single event has a much greater probability of being back- 
ground than at 10 TeV, where the source signal is dominant. 
Spectra allow for such an approach, which takes full advan- 
tage of experimental data, in studying the high energy behav- 
ior of TeV sources. Resolving hadronic activity at these ex- 
treme energies would provide clear evidence concerning the 
sources of Galactic cosmic rays lll3ll . 

Atmospheric cosmic-ray showers give rise to high rates of 
down-going muons, forcing a neutrino telescope to search for 
up-going muons resulting from neutrino interactions. Ice- 
Cube is well-situated to utilize the high resolution of these V^- 
induced muons in observing northern-sky sources. However, 
a detector is needed in the northern hemisphere to accurately 
locate southern-sky neutrino sources, although IceCube may 
also measure shower rates from particularly bright sources. 
Together, IceCube and a km 3 Mediterranean detector will pro- 
vide full-time coverage of the entire sky, a feature distinct to 
neutrino telescopes. The combined observations from these 
detectors can be used to study compound objects, like the Vela 
complex, by: (1) Discovering neutrino sources through high- 
resolution Vu-induced muons; (2) Confirming agreement with 
gamma-ray observations in the low energy regime; (3) Exam- 
ining previously unexplored energies using muons and show- 
ers together. 

The sources that we examine in detail are briefly discussed 
in Section|n] Our empirical calculations, based upon observed 
TeV spectra, should be compared with previous Galactic neu- 
trino studies, which have considered theoretical production 
models or less recent data (e.g. Refs Q 13 111 O El E 
|H Ell |13 S). In two very recent works, byLi- 
pari 1 8] and Kappes, Hinton, Stegmann, and Aharonian |9], a 
similar approach to ours has been followed, with complemen- 
tary perspectives and details. 



H. PROMISING TEV NEUTRINO SOURCES 

Vela Region — Of particular interest amongst prospective 
neutrino sources, the shell-type supernova remnant (SNR) 
Vela Jr. (RX J0852.0-4622) is one of the brightest objects 
in the southern TeV sky. The hard, intense TeV gamma-ray 
spectrum of Vela Jr., best explained as being pionic in nature, 
makes it an intriguing object to study |26]. Additionally, the 
morphology of the TeV emission may make this one of the 
most interesting astrophysical neutrino sources. HESS has 
also measured a TeV spectrum from Vela X, the pulsar wind 
nebula (PWN) associated with the larger Vela SNR, and ad- 
vanced a leptonic origin l27ll . However, if this spectrum is 
instead pionic, as proposed in Ref. 12811 . the accompanying 
neutrino flux would be easily detectable. The non-detection 
of neutrino events from such an intense source would allow 
for a significant test of leptonic production in a short period of 
observation. 

Galactic Center Diffuse Emission — HESS has also re- 
cently discovered a region of diffuse TeV emission from the 
Galactic Center ridge |29|. This source has several interesting 
features which merit further investigation. The large extent 
of the emission, hardness of the spectrum, high gas density 
(which is well-correlated with the TeV emission) and strong 
magnetic fields in the region leave very little doubt that this 
spectrum is pionic l29Tl . In addition, the total flux from this 
region is actually about twice as intense as that of the pre- 
viously discovered source coincident with Sgr A* [30]. Mea- 
surement of the accompanying neutrino flux would provide an 
independent confirmation of the means of production, some- 
thing which has never been possible. 

Unknown Knowns and Known Unknowns — We also ex- 
amine the neutrino detection prospects of other known TeV 
sources, including SNR RX J 17 13-3946, the spectrum of 
which has been measured twice by HESS and determined to 
likely be pionic 13 ll 13211 . The calculated Vu -induced muon 
fluxes based upon the two different sets of HESS data il- 
lustrate the effects of different high energy spectral assump- 
tions. Recently, HESS has reported the discovery of four TeV 
sources in the Galactic Plane which have no apparent coun- 
terparts at other wavelengths 13 3l 13411 . In addition to these, 
EGRET previously discovered numerous unidentified sources 
of GeV gamma rays in the Milky Way 1 35 ] . The confirmed 
observation of neutrinos from these sources would also reveal 
an intimate connection between TeV gamma rays and cos- 
mic rays. Furthermore, the understanding of these and other 
gamma-ray sources may account for the TeV excess |36]. 



m. BASICS OF NEUTRINO DETECTION 

In high energy p-p scattering, 7T + , 7t~, and 71° are pro- 
duced in nearly equal numbers 13711 . Gamma rays are the 
result of the decay K° — > 77, while neutrinos originate from 
the 7T+ -> jU+V M -> e+v^VeV^ and vT -> jll^ -> e~V M V e V M 
decay channels. The ratio of neutrinos to photons from pion 
decay is easily found. From charged pion decay, the resulting 
initial neutrino flavor ratio, v e : : V T , is 1 : 2 : 0. During the 
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traversal of astrophysical distances, vacuum neutrino oscilla- 
tions transform this ratio to 1 : 1 : 1 . In neutrino telescopes, 
neutrinos and antineutrinos are practically indistinguishable. 
We can then consider their sum, V + V, and average cross sec- 
tion in all calculations. All further references to neutrinos will 
imply v + V. Thus, for equal pion multiplicities, each pho- 
ton from 7T° decay corresponds to one neutrino of each fla- 
vor (N r = N Ve = N VlA = N Vt =N v ). The typical energy of the 
neutrinos resulting from these decays is ~ 1/2 of the gamma- 
ray energy from n° decay. The resulting v + V spectrum is 
then shifted, relative to the source gamma-ray spectrum of 
dQyjdEy = <j) r E~ r , as 



d<$> v 

dE v 



r-i 



(j) 7 E v r = <p v E 



-r 



(1) 



where we consider each neutrino flavor separately. 

We will generally refer to the normalization of the gamma- 
ray spectrum as <py, the differential photon flux at 1 TeV (in 
TeV~' cm~ 2 s _1 ), which is commonly used in gamma-ray 
astronomy. Typical source gamma-ray spectra have (j)y ~ 5 — 
20 x 10~ 12 TeV- 1 cm" 2 s" 1 and F ~ 1.8-2.4. For greater 
detail on the relationship between the initial proton spectrum 
and the resulting product spectra, see Ref. 1 38 ] . Pions can 
also be produced, in a different ratio, in p-y scattering, but 
this is generally only important at much higher energies |39]. 
A neutrino flux may be observed through neutrino-induced 
muons and electromagnetic (and hadronic) cascades, which 
are referred to as showers. We shall consider both methods in 
further detail. 



A. Muon Detection 

Our analysis of muon detection will be limited to Vj, 
charged-current (CC) events with an observable final en- 
ergy of Efx > 0.5 TeV, which can be detected through their 
Cerenkov emission with an angular resolution of < 1° i40ll . 
It is important to recognize that the relevant quantity in these 
events is the measured energy of the muon, which can be re- 
constructed from radiative losses in the detector l4ll l42l l43l 
l44ll . As muons can be produced far outside of the detector, 
with energy loss prior to entering, results of calculations are 
given in terms of E^ (at the detector) instead of E v . 

We can consider -induced muons produced both inside 
and outside of the instrumented volume of the detector. Con- 
tained event rates can be found by combining the neutrino 
flux, detector mass, and afcc f° r neutrino-nucleon scattering. 
We use the average of the v-N and v-N cross sections, as 
computed in Ref. |45], and (y(E v )) from Ref. |46]. We con- 
sider a detector that is entirely composed of, and enclosed by, 
water (ice), with a km 2 detector area, which is a reasonable 
approximation at these energies. Using an IceCube-like ef- 
fective detector area, such as computed in Ref. [47], would 
lead to the lower energy muon spectrum being slightly sup- 
pressed. An actual km 3 detector will be situated on top of 
solid rock, which will enhance the upward-going neutrino- 
induced muon rates relative to our calculations. Also, muons 



produced via v T CC interactions, which contribute to the to- 
tal muon flux through the T — > /i ± V fl V T decay channel, will 
not be considered 1371 l48l l49ll . Our formulas follow those 
of Gaisser l37l Boll , with appropriate approximations. Con- 
sidering the uncertainties in both the gamma-ray spectra and 
neutrino telescope performance, which may not be fully un- 
derstood until the actual detectors are calibrated, as well as 
the low statistics, our calculations are sufficient for the scope 
of this work. 

Using these parameters, the spectrum of contained V^- 
induced muon events is calculated as 



dE u 



KV A J-^e- E ^ acc {E v )e-^. 
dE v 



(2) 



where Vd et is the detector volume, Ey is an assumed ex- 
ponential cutoff in the neutrino spectrum, and the term K = 
NaPT (1 — y(E v )) ~ ! takes into account observation time (T), 
normalization of the muon spectrum, and the molar density of 
water. The energy of the produced muon is related to the orig- 
inal neutrino energy as = (1 — y(E v ))E v . The term e~ T *, 
with T e = Na^b °tot(Ev)> accounts for neutrino attenuation 
due to scattering within Earth. Here, o tot = dec + Onc and 
A® is the average column depth (in cm.w.e.) of Earth based 
upon the declination of the source 1 46] . This attenuation fac- 
tor only becomes important at Ey> 10 TeV, and varies with 
declination, as discussed in Ref. |8]. 

A high energy muon born outside of the detector may still 
be detectable when it enters the instrumented detector volume. 
To find the propagation range of a muon, we assume an aver- 
age continuous muon energy loss of 



dE 
dX 



-a-fiE, 



(3) 



where a = 2.0 x 10~ 6 TeV cm 2 g" 1 and j3 = 4.2 x 10~ 6 cm 2 
g _1 l42l l48ll . Integrating the average energy loss results in a 
range for a muon of initial energy E^ of 



1 



R^E^Ef,) = -In 



a + pEl 



(4) 



where E^ is the energy of the muon as it enters the detector. 
As the muon range (typically a few km) increases with energy, 
so does the effective volume of the detector. The observed 
through-going spectrum, accounting for the probability of ob- 
serving a muon entering the detector with energy En, can then 
be calculated as 



dN„ 
dE u 



thru 



N A pTA det 



a 



E u aE v 



e - £ v/£$ Ut acc (£ v ) e -% )(5) 



where A^et is assumed to be 1 km 2 . It should be noted that 
our use of the total afcc and (y) may overestimate the muon 
signal rates by ~ 20 — 30% (perhaps more for spectra harder 
than E~ 2 ). However, the atmospheric background will also be 
overestimated by a similar, even slightly larger amount. 
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B. Shower Detection 

For bright point sources or low surface brightness extended 
sources, shower events caused by V e , V T CC interactions can 
be used to effectively determine the source spectrum. To de- 
tect shower events, they must be at least partially enclosed 
within the instrumented volume. The shower signal is dom- 
inated by v e and V T CC events, which transfer ~ 100% of 
the original neutrino energy to the shower (E^ « E v ). The 
measurable shower spectrum, given in terms of the observed 
shower energy, E s \„ 

=2N AP TV d J-^e- E ^o cc (E v )e-^, (6) 



\dE sh J 



CC 



dE 



then effectively traces the original neutrino flux, as the 
neutral-current (NC) events from all three flavors, which only 
contain E s h = (yNc)E v (in addition to the smaller Onc), only 
account for < 10% of the total signal. 

Combined shower and muon observations would deter- 
mine the V^/(v e + V T ) ratio arriving from a source. This 
measured flavor ratio has many applications, as detailed in 
Refs. (i3 EH EH |53l . The atmospheric v e flux is only about 
~ 1/20 of the Vji background, which itself only adds to the 
shower background through the weaker NC channel, since 
CC interactions are identifiable by the resulting muon tracks. 
We will not consider the prompt background, which is the 
only source of atmospheric V T , as it only becomes important 
at very high energies 15 311 . These facts help offset the lower 
angular resolution of these events (< 10° is ex pect ed for a 
Mediterranean detector |54], < 25° for IceCube |40]), which 
sets the background event rate. The improved shower reso- 
lution of a Mediterranean detector, as compared to IceCube 
(which is naturally limited by light scattering in ice), increases 
the value of showers in determining the high energy properties 
of a source. 



C. Neutrino Spectroscopy 

To demonstrate the ability of a neutrino telescope to differ- 
entiate between various spectral properties, we first consider 
a general object with a fixed <j)y = 20 x 10~ 12 TeV 1 cm~ 2 
s _1 (at 1 TeV) and assume a source declination of 8 = +10°. 
The resulting up-going muon fluxes for a range of spectral 
indices and high energy exponential cutoffs, for one year of 
IceCube observation, are shown in Fig. [2] In the top panel, 
we consider spectral indices with r= 1.8, 2.0, 2.2, and 2.4, 
with Ey Ut = 50 TeV. Note that the differences are not entirely 
due to differing integrated neutrino fluxes, as the rising <7 V 
results in higher event rates for harder spectra. In Ref. 
similar results were obtained by normalizing to the integrated 
flux above 1 TeV. The bottom panel illustrates the effect of 
100, 50, 25, and 10 TeV neutrino spectrum exponential cut- 
offs for a fixed index (r = 2.1). IceCube is well-suited to iso- 
late a similar point source located in the northern sky through 

-induced muons within a short period of operation. Such 
a source would also be observable through shower events by 



a km 3 Mediterranean detector. The improved shower detec- 
tion capabilities of a Mediterranean detector would allow for 
a more direct measurement of the arriving neutrino spectrum, 
adding greatly to the science that can be extracted from TeV 
sources. 

All of our calculations are done on an empirical basis us- 
ing the average measured parameters of gamma-ray source 
spectra. All high energy exponential cutoffs that we will 
consider are given in terms of Ey Ut (=l/2Ey Ut ), to allow for 
a range of calculations to be compared with future experi- 
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FIG. 2: A demonstration of the sensitivity of neutrino telescopes to 
varying spectral characteristics. Shown is the (base- 10 log) differ- 
ential (Vjj + v^)-induced muon spectra for a pionic source spectrum 
with (fly = 20 x 10~ 12 TeV -1 cm~ 2 s~' . In the upper panel, the lines 
of decreasing height correspond to spectral indices of T = 1.8, 2.0, 
2.2, and 2.4, with an exponential cutoff Zs£ ut = 50 TeV. In the lower 
panel, the lines of decreasing height correspond to neutrino spec- 
trum cutoffs of 100, 50, 25, and 10 TeV, with a fixed spectral index 
of r = 2. 1 . Each line corresponds to one year of observation. 
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mental data. Angle-averaged FLUKA computations of at- 
mospheric neutrinos past E v = 1 TeV are used to determine 
backgrounds ll55ll . For ease of comparing the expected sig- 
nal to background and assessing total rates, we will mainly 
provide integrated rates above a given measured energy. Be- 
cause the atmospheric spectrum is very steeply falling at these 
energies (~ £y 3,s ), it is often advantageous to use a low en- 
ergy event cutoff > 1 TeV when calculating significance. Our 
method can also yield the measured muon and shower spectra 
directly. Spectral information should be used for maximum 
likelihood analysis of observed neutrino events in a neutrino 
telescope, to make the best use of measured energies in deter- 
mining detection significance. This spectral information re- 
quires consideration of the energy resolution of the detector. 
In IceCube, the energy resolution for muons and showers is 
expected to be ~ 20% and ~ 10% in the logarithm of the en- 
ergy, respectively 1 40, 56], which does not significantly affect 
our conclusions. 

Note that we present our results for an assumed muon 
effective area of 1 km 2 , and assumed muon angular reso- 
lution of < 1°, both independent of energy. These results 
are close to the results of detector simulations, especially > 
1 TeV, the primary area of interest 1 40, 57]. Below 1 TeV, 
these assumptions are certainly overly optimistic; however, 
in that range, the atmospheric background is dominant, and 
so accuracy is of less importance. In regards to the an- 
gular resolution, our assumptions are somewhat too conser- 
vative (in particular, the improved resolution at higher en- 
ergy will help reject background). Part of our intent is to 
show the likely signals, and to encourage the experimental- 
ists to optimize their detector designs in order to achieve 
the necessary sensitivity to observe Galactic TeV sources. 
These calculations are not meant to replace a more detailed 
study of neutrino source detectability, which would include 
time-dependent source locations, zenith-angle-dependent at- 
mospheric backgrounds, angle-dependent detector sensitivi- 
ties, event reconstruction methodologies, stochastic muon en- 
ergy losses, source modeling, etc. Each prospective source 
should be subjected to such a full Monte Carlo simulation by 
the IceCube and Mediterranean collaborations. 



IV. THE VELA COMPLEX 

Of the many known TeV gamma-ray sources, the shell- 
type SNR Vela Jr. (RX J0852.0-4622) is one of the most 
interesting. This southern-sky source has been observed 
at gamma-ray energies exceeding 10 TeV by HESS. Their 
analysis is suggestive of a hadronic origin for the gamma- 
ray spectrum 12611 . This source is very bright in gamma 
rays {d<S>/dE = 21 x Kr^/TeV) -2 - 1 TeV" 1 cm" 2 s" 1 ), 
with well-defined regions of gamma-ray emission. Shell-type 
SNRs are considered to be the most likely sites of Galac- 
tic cosmic -ray proton acceleration |58l]. As these source 
proton spectra are expected to be cut off near the knee (~ 
3 x 10 15 eV), a cutoff should also be present for gamma-rays 
and neutrinos at a lower energy scale 13811 . Thus, we calculate 
the expected neutrino-induced muon rate assuming a pionic 



spectrum with several neutrino spectrum exponential cutoffs 
(50, 25, and 10 TeV), as shown in Fig. [2] When considering 
emission from the entire source extension, we must accept at- 
mospheric background from a ~ 7 deg 2 area. Comparing the 
source and background muon rates, this SNR is expected to 
be significantly detectable with a km 3 Mediterranean detec- 
tor. For the higher cutoffs, discovery may be possible in only 
a few years. These event rates also suggest that ANTARES 
may be able to find evidence of neutrino emission from this 
source |59]. 

A number of features make Vela Jr. unique among prospec- 
tive TeV neutrino sources. The TeV emission is observed to 
originate from several regions which are separated by ~ 2° . In 
addition to detecting the presence of neutrinos, the intensity of 
the measured muon spectrum may provide adequate statistics 
to allow accurate location of neutrino-production sites. Thus, 
this unique gamma-ray morphology allows for the possible 
construction of a neutrino map of Vela Jr. This capability 
would enable a comparison of the production mechanism and 
high energy activity in different regions of the same source. 
When observing features of the neutrino source, the total at- 
mospheric neutrino background enclosed is lower, which in- 
creases their individual detectability. An example of such a 
map can be seen in Ref. |60]. Further HESS observations 
should provide improved resolution of the TeV emission. 

Additionally, the brightness of the source may also make 
it possible to observe appreciable numbers of showers. As a 
consequence of the detected showers directly tracing the neu- 
trino flux, spectral information may be found (even with a 
slightly higher background). Fig.|4]shows the rate of showers 
for one year of observation, compared to the irreducible atmo- 
spheric background. With enough observation time, it may be 
possible to determine the neutrino flavor ratio from this source 
by observing the spectrum with a km 3 Mediterranean de- 
tector and by also utilizing IceCube to effectively increase the 
shower volume (since shower events must be located within 
the detector) in measuring the v e + V T flux. 

Vela Jr. is coincident with a region of the larger Vela SNR. 
HESS has also reported the discovery of TeV gamma rays 
originating from the Vela PWN, referred to as Vela X (of no 
relation to the HMXB Vela X-l) ||3|. Vela X has a very hard 
spectrum, which is best fit with a spectral index of T = 1 .45 
with an exponential cutoff beginning at Ey Ut = 13.8 TeV. In 
their discovery paper, HESS concludes that this emission is 
produced by an inverse Compton mechanism. An independent 
analysis has considered the possibility that this spectrum is 
produced hadronically, with an accompanying flux of neutri- 
nos lE8ll . This source is sufficiently separated from Vela Jr. to 
independently test for neutrino emission through muons. If 
we assume the entire measured flux is pionic, then the ex- 
pected -induced muon rate is N(E^ > 1 TeV) ~ 4.5 yr _1 , 
in a km 3 Mediterranean detector, which would be relatively 
easy to measure. While our calculation may somewhat over- 
state the muon rate for such a hard spectrum, nevertheless, 
the nondetection of a significant neutrino event excess would 
indicate leptonic production. 

Considering the large apparent size of the Vela SNR shell 
(~ 8°), any TeV gamma-ray emission originating from it may 
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FIG. 3: Vela Jr. (RX J0852.0-4622)—Muons: Integrated (v M + 
\>H )-induced muon rates above a given measured muon energy. The 
long dashed (blue), dot-dashed (red), and short dashed (gold) lines 
correspond to neutrino spectrum exponential cutoffs of 50, 25, and 
10 TeV, respectively. The dotted line shows the expected atmospheric 
background in a 7 deg 2 bin, as discussed in the text. Rates are for one 
year of operation in a km 3 Mediterranean detector. 



be too diffuse to detect directly. Searching for it in neutrinos 
is a possibility. A northern-hemisphere detector would isolate 
any neutrino point source spectra located in the region through 
V^-induced muons. Showers have the ability to observe this 
entire region simultaneously, raising the intriguing prospect of 
utilizing IceCube and a km 3 Mediterranean detector in concert 
to discover and study bright, but highly extended TeV neutrino 
sources. 



V. GALACTIC CENTER REGION 

The HESS discovery 12911 of the region of Galactic Cen- 
ter diffuse emission (GCD) is important for neutrino, gamma- 
ray, and cosmic -ray astrophysics 16 U 16211 . The spectrum of 
the GCD was measured over a very large region spanning 
the Galactic coordinates |/J < 0.8°, \b\ < 0.3° with d<$/dE = 
1.73 x 1(T 8 (£/TeV)- 2 - 29 TeV" 1 cirr 2 s" 1 sr" 1 . A num- 
ber of large, dense (hh ~ 10 4 cm~ 3 ) molecular clouds are 
known to fill this region |63]. Considering the vast extent 
of this hard emission (~ few hundred pc across), along with 
the high target density and magnetic fields in the region, the 
only reasonable production mechanism for these gamma rays 
is neutral pion decay. The close morphological correlation 
between the TeV emission and the gas distribution is particu- 
larly compelling 1 29] . This TeV gamma-ray spectrum implies 
a local cosmic-ray spectrum that is much harder and 3 — 9 
times denser than that of the CR flux as measured at Earth. 



FIG. 4: Vela Jr. (RX J '0852. 0-4622)— Showers: Integrated 
neutrino-induced shower rates above a given measured shower en- 
ergy, with signal lines as in Fig. [5] Note that a different energy scale 
is used for showers, as measurable shower events are assumed to 
have a 1 TeV energy threshold. Background rates are for a circle 
of 10° radius, corresponding to the shower angular resolution in a 
km 3 Mediterranean detector. A significant reduction in background 
results from improved angular resolution. Rates are for one year of 
operation. 



The likely source of these cosmic-ray protons is Sgr A, which 
hosts the remnant of a recent supernova Ir34 |l63l . These pio- 
neering conclusions were reached in Ref. 12911 . 

In order to calculate the expected -induced muon rate, we 
must first take into account the extent of the GCD source. The 
muon angular resolution of a neutrino telescope at these en- 
ergies (~ 1°) covers the entirety of the GCD. In effect, this 
region can be treated as a neutrino point source. Integrat- 
ing over the ~ 1 deg 2 source region, the total photon flux is 
d<P/dE = 5.2 x Kr 12 (£/TeV)- 2 - 29 TeV" 1 cirr 2 s" 1 . We 
define the differential flux at 1 TeV as 0gcd = 5.2 x 10~ 12 
TeV -1 cm~ 2 s~'. We will not consider showers from the 
GCD; however, any such study should also take into account 
a similar diffuse emission around Sgr B 1 29 ] . 

We first consider the conservative assumption that the GCD 
spectrum exhibits an exponential neutrino spectrum cutoff 
Ey Ut = 20 TeV. As a cosmic -ray production site, we also con- 
sider a higher cutoff of 50 TeV. The expected muon rates 
are shown in Fig. [5] Even though its angle-integrated flux 
is larger, the GCD was revealed only after subtracting the 
previously discovered TeV emission from Sgr A* and a dis- 
tant SNR HO, EH. The TeV gamma-ray flux from Sgr A*, 
which has <py ~ 0.5 <Pgcd, can reasonably be assumed to also 
be pionic, as the spectral index (r = 2.2) is close to the GCD 
value 1 30, 67]. As the two sources are coincident, we combine 
their expected muon rates, as also shown in Fig. [5] Such an 
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approach to calculating the total neutrino flux was also taken 
in Ref. 1 68], where the integrated muon rate was found. The 
declination of the source, along with the latitude of the de- 
tector, determine the fraction of the time that the source is 
below the horizon 12411 . For the Galactic Center this is ~ 0.7, 
which is taken into account in Fig. [5] The signal rates shown 
in Fig. [5] when compared to the expected number of atmo- 
spheric events, suggest that this source is near the threshold 
of discovery for the lifetime of a km 3 Mediterranean detector. 
Discovery may be hastened by utilizing maximum likelihood 
methods in comparing the expected signal with background 
for each muon event with a measured energy, as higher energy 
events have a much greater significance. Such an analysis is 
only possible with spectral information. 



VI. A MULTITUDE OF TEV SOURCES 

Numerous discoveries of new gamma-ray sources should 
be expected from the HESS, MAGIC and VERITAS air 
Cerenkov telescopes and the GLAST satellite l69l FTol l7lll . 
However, the characteristics of the source spectra above ~ 
10 TeV are difficult to observe. RX J1713-3946 is also a 
shell-type SNR, first observed in the TeV regime by CAN- 
GAROO 0]. It has since been observed twice by HESS 
with greatly improved statistics. We consider both HESS 
observations to illustrate the potential of neutrino telescopes 
to determine the high energy behavior of TeV sources. The 
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FIG. 5: Galactic Center region: Integrated (Vh + v^)-induced muon 
rates above a given measured muon energy. Shown are the rates from 
the Galactic Center diffuse source (GCD) alone (dot-dashed) and the 
result of including the Sgr A* source (dashed). Two cases are con- 
sidered: a neutrino spectrum exponential cutoff £y Ut = 50 TeV and 
a more conservative cutoff of 20 TeV. The dotted line shows the es- 
timated atmospheric background in a 3 deg 2 bin. Rates are for one 
year of observation in a km 3 Mediterranean detector. 



flux was initially reported up to 10 TeV as a power law 
(T ~ 2.19), with analysis favoring hadronic production 1 3 1 ] . 
Subsequent observations revealed a spectrum that extends to 
at least 40 TeV with steepening past ~ 10 TeV that can be fit 
by either an energy dependent spectral index or an exponen- 
tial cutoff 113211 . Previous neutrino studies of this source were 
inspired first by the CANGAROO detection and later by 
utilizing the more accurate 2004 HESS spectrum along with 
a more detailed calculation fl8ll . We analyze the -induced 
muon spectra from these two data sets independently, assum- 
ing E£ ut = 50 TeV for the 2004 data. As seen in TableUH the 
integrated rates above 1 TeV are similar. However, at higher 
energies calculations based upon the old power law spectrum 
exhibit a deviation from the new data. From the 2004 data, we 
find 7V(£ M > 10 TeV) - 0.3 yr"\ while the 2005 data yields 
a rate only ~ 1/10 of this, which highlights the sensitivity of 
neutrino observations to the spectrum at the highest energies. 
Based on either spectrum, this SNR should be detectable in a 
km 3 Mediterranean detector. 

Amongst the HESS catalog of TeV sources are several with 
no observed counterpart at any other wavelength. Explaining 
the origin of these sources is a theoretical challenge. It is pos- 
sible that they represent a non-trivial new class of astrophysi- 
cal objects. With so little being known about them at present, 
the observation of neutrinos from these sources would yield 
invaluable insight into their nature. The first of these sources 
is HESS J1303-631 [34]. Chandra observations revealed no 
likely x-ray counterpart |73]. This source is always visible to 
a Mediterranean detector. Three sources were also found in 
the HESS survey of the Galactic Plane that remain unidenti- 
fied: HESS J1614-518, HESS J1702-420, and HESS J1708- 
410 m. The spectral information and resulting integrated 
V)i -induced muon rates from these sources are given in Ta- 
ble|I] It has been proposed that HESS J 1303-631 is a remnant 
of a Galactic gamma-ray burst (GRB), with a harder (r ~ 2.2) 
and more extended emission than has been observed Q. 
From the given predictions of the proposed hadronic source 
model, we would expect N(E^ > 1 TeV) ~ 2.1 yr , which 
is potentially verifiable with a km 3 Mediterranean detector. 
From the measured spectra, detecting these sources individu- 
ally may be difficult; however, the stacking method may prove 
advantageous. Stacking effectively increases the total obser- 
vation rate by combining the individual event rates from iden- 
tical sources 11751. 17611771 17811 . Stacking of the observed muon 
data may improve the prospects of determining whether neu- 
trinos are present and help unveil the mystery of these un- 
known sources. 

Other Galactic TeV sources with probable counterparts may 
also be examined as potential neutrino sources. Many new 
sources were found as a result of the HESS inner Galactic 
Plane survey, some of which may be SNRs 13311 . MAGIC 
and VERITAS observations will reveal sources in the outer 
regions of the Plane, which has yet to be examined with an 
instrument of their sensitivity. Assuming that the gamma-ray 
spectra already observed by HESS are pionic, the calculated 
neutrino event rates for a number of sources are shown in Ta- 
ble mi Several of these have sufficient rates to consider as 
possible neutrino sources independently. However, it would 
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TABLE I: Unidentified HESS Sources: Integrated (v^ + Vfi)-induced 
muon rates, assuming a pionic spectrum, for a 1 TeV threshold 
and neutrino spectrum exponential cutoff E"". For one year of km 3 
Mediterranean detector operation, accounting for observable time 
below the horizon. The HESS differential flux at 1 TeV is given 
in terms of 10~ 12 TeV -1 cm~ 2 s _1 , with spectral index F. 



Source 


<t>7 


r 


E c v l " (TeV) 


AV(> ITeV) 


HESS J1303-631 


4.3 


2.44 


10 


0.3 








50 


0.5 


HESS J1614-518 


7.0 


2.46 


10 


0.5 








50 


0.8 


HESS J1702-420 


2.5 


2.31 


25 


0.3 








50 


0.4 


HESS J1708-410 


1.5 


2.34 


10 


0.1 








50 


0.2 



be difficult to significantly detect dimmer sources over the at- 
mospheric background. Stacking, in combination with a max- 
imum likelihood spectral analysis, may be used to examine 
the entire group, or subgroups, of the lower-flux sources. 

VII. A TEV NEUTRINO BEACON 

It would be beneficial to have a bright source of TeV neu- 
trinos that can be seen by both IceCube and a Mediterranean 
detector. The ideal source for this purpose would have a decli- 
nation in the range + 10° < 8 < +30°, which also has the ben- 
efit of reduced neutrino attenuation in Earth at high energies, 
resulting in increased muon and shower event rates. The Crab 
nebula is already an accepted standard candle in TeV gamma- 
ray astronomy. From HESS measurements 1 79], if the entire 
TeV spectrum could be attributed to hadronic processes, our 
methods yield an event rate of N(E fl > 1 TeV) ~ 2.7 yr~' in 
IceCube. The Crab spectrum is relatively soft and well de- 
scribed by leptonic processes [80], although AMANDA ob- 
servations hint at a neutrino signal llsill . A pionic compo- 
nent may be uncovered by a significant observation of neu- 
trino events. 

The Vela complex that we have already discussed is located 
in the southern sky. Only after extensive observations was the 
existence of Vela Jr. confirmed |82]. It would not be surpris- 
ing, then, for a similarly bright TeV source to be discovered 
in the northern sky (perhaps directly through neutrino events). 
The spectrum for such an object was calculated near the end 
of Section[ni]for a variety of spectral indices and cutoffs. Fu- 
ture gamma-ray observations of this declination range in the 
TeV regime may reveal another important source for neutrino 
astronomy. 

Multi-wavelength observations of such a source would add 
a rung to a TeV neutrino "distance ladder". Cosmic -ray inter- 
actions in the Earth's atmosphere generate the nearest guar- 
anteed neutrino source 1831 18411. Similar processes in the at- 
mosphere of the Sun, which is at a well known distance, are 
quite likely to produce TeV neutrinos as well IsEl. l86l l87ll . 



TABLE II: HESS Sources With Probable Counterparts: Integrated 
(Vfx + vy,)-induced muon rates, assuming a pionic spectrum, for a 
1 TeV E^ threshold and neutrino spectrum exponential cutoff E"" . 
For one year of km 3 Mediterranean detector operation (unless noted), 
accounting for observable time below the horizon. The HESS differ- 
ential flux at 1 TeV is given in terms of 10~ 12 TeV~ 1 cm~ 2 s~ 1 , with 
spectral index F. 
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Of the prospective sources, the Vela SNR is estimated to be 
only ~ 250 — 300 pc from Earth, while Vela Jr. may be even 
closer i8| j89ll9(il . RX J1713-3946 has a distance estimate of 
~ 1 kpc 19111 . At ~ 8 kpc, the GCD is the most distant guaran- 
teed source with a confirmed distance. AGN, if they produce 
measurable fluxes of neutrinos, would provide a range of very 
remote point sources. Measurements of neutrinos from a vari- 
ety of distances (ideally with flavor ratios) would provide im- 
portant information for testing neutrino properties at increas- 
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ing L /E (distance/energy) 
100]. 



VIII. EGRET, THE TEV EXCESS, AND THE GREAT 
UNKNOWN 



The third EGRET catalog contains many uniden tifie d 
gamma-ray sources detected in the MeV-GeV range 

Gnu. 

In the Whipple survey of a selected group of such EGRET 
sources, upper limits on high energy gamm a-ray emissions 
from a number of sources were determined ill 0211 . Two of 
these sources, 3 EG J 1337+5029 and 3EG J2227+6122, pos- 
sessed measured excesses that were suggestive of gamma-ray 
emission 1 103]. The upper limits placed on these sources do 
not significantly constrain the EGRET data. We consider the 
possibility that these spectra are pionic and continue into the 
TeV regime. Spectral information, cutoffs, and muon rates 
are given in Table [II]] These would be readily detectable 
by IceCube. We note that 3EG J2227+6122 is roughly co- 
incident with a position in the three-year AMANDA point 
source survey which exhibits a modest excess of neutrino 
events [12]. A less significant excess still appears in newer 
data 1 104]. In general, a maximum likelihood analysis assum- 
ing the spectrum and normalization of such bright EGRET 
sources extended into the TeV range may prove useful for neu- 
trino searches. Future observations by GLAST, VERITAS, 
and MAGIC will provide valuable additional spectral infor- 
mation. 

Recently, Milagro announced the detection of diffuse TeV 
emission from a ~ 600 deg 2 patch of the Galactic Plane that 
contains ten unidentified EGRET point sources |105|. This 
measurement places high energy constraints on the spectra of 
these ten sources, as discussed in Ref. |36]. Consequently, 
these spectra, many of which are quite hard, cannot all simply 
extrapolate to TeV energies. HEGRA discovered an uniden- 
tified TeV source in this region, however, it ca n account for 
only a small fraction of the Milagro flux 110611 . Located in 
the northern sky, this region is accessible to observations by 
IceCube. If this emission is pionic and truly diffuse, then 
it woul d be very difficult to detect through neutrino-induced 
muons 1 107]. Assuming that the entire flux is instead divided 
among several point sources, the expected muon rate can be 
calculated. A detailed Milagro study of the Cygnus region 
revealed a significant TeV flux 1 108]. Assigning all of the Mi- 
lagro 3.5 TeV flux to a single source and extrapolating with a 
typical r — 2.2 gamma-ray spectrum (out to 100 TeV), we ar- 
rive at a limiting case for the possible neutrino emission from 
the region. For a pionic spectrum, the calculated -induced 
muon rate is N(E )1 > 1 TeV) ^2.3 yr~ ! in IceCube. Fu- 
ture measurements of this region by Milagro, GLAST, VER- 
ITAS, and MAGIC will determine whether multiple sources 
contribute to this diffuse flux. 



TABLE III: Two Unidentified EGRET Sources: Integrated (v p + V M )~ 
induced muon rates, assuming a pionic spectrum, for a given E^ 
threshold and neutrino spectrum exponential cutoff E"" (in TeV). 
For one year of IceCube operation. The assumed differential pho- 
ton flux at 1 TeV is given in terms of 10~ 12 TeV -1 cm~ 2 s _1 , with 
spectral index T. 



Source 



r gg 



0.5 TeV) jV M (>lTeV) 



3EG J 1337+5029 20 2.2 



3EG J2227+6122 10 2.0 



10 
50 
10 
50 



4.1 

6.9 
3.0 
5.5 



2.5 
4.7 
1.9 
4.0 



IX. DISCUSSION AND CONCLUSIONS 

Upcoming neutrino telescopes will deliver the first di- 
rect evidence concerning the production mechanism of 

TeV gamma-ray sources. Following the pioneering 

AMANDA lll09ll and Baikal ill ldL efforts, the next generation 
IceCube |6] and Mediterranean |7] km 3 detectors will reach 
the scale necessary to examine the Galactic sources discussed 
here. Hadronic mechanisms would be confirmed through neu- 
trino detection, while a significant absence of neutrino events 
would imply leptonic processes. They will be able to probe 
the expected high energy spectral cutoffs that would other- 
wise be unobservable. It has been suggested that the cutoff in 
a SNR may evolve with age | 111]. Finding these cutoffs with 
neutrinos may also yield valuable new information in this re- 
gard. 

Combined, IceCube and a km 3 Mediterranean detector will 
provide continuous, all-sky coverage and can be used together 
to detect neutrino-induced muons and showers from TeV 
sources. The good angular resolution for muon events can be 
used to precisely locate a source in neutrinos. Showers, and 
contained muon events, provide accurate reconstruction of the 
source spectrum at energies beyond the reach of gamma-ray 
telescopes. These measurements will provide important infor- 
mation concerning the origin of high energy Galactic cosmic 
rays, potentially directly observing the source population re- 
sponsible for production up to the knee at ~ 3 x 10 15 eV. 

The shell-type SNR Vela Jr. (RX J0852.0^622) is an 
intriguing prospective source of TeV neutrinos. The ex- 
tent and intensity of the TeV emission make this, as well 
as the Vela complex as a whole, a unique target for neu- 
trino telescopes. The confirmation of hadronic/leptonic pro- 
cesses through detection/non-detection of neutrino events 
from Vela X and the possible detection of neutrinos through 
shower events from the Vela SNR shell are exciting possibili- 
ties. 

While it is expected that other sources may eventually be 
more compelling, only the Galactic Center diffuse emission 
can presently be claimed to possess a guaranteed neutrino 
flux. The detection of this neutrino flux will give further 
insight into the complex processes occurring in the Galactic 
Center region. When we consider that, to date, no high energy 
astrophysical neutrinos have ever been positively detected, the 
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importance of such a guaranteed flux of TeV neutrinos is dif- 
ficult to overstate. 

Included in the catalog of TeV sources are several that 
remain unidentified. As more observations are undertaken 
by the HESS, MAGIC, and VERITAS telescopes, this num- 
ber is expected to increase. Discovering neutrino fluxes 
from these sources would provide invaluable information con- 
cerning their nature. As new TeV gamma-ray sources are 
discovered, the pool of potential TeV neutrino sources in- 
creases. While some of these may not be significant alone, 
when grouped into classes, stacking may increase the poten- 
tial for neutrino studies. An even larger number of unidenti- 
fied sources remains in the GeV regime from EGRET obser- 
vations. Many of these sources have intense, hard spectra. If 
any of these sources, some of which are visible to IceCube, 
have spectra that extrapolate into the TeV regime and are pio- 
nic, they would have abundant fluxes of neutrinos. 

In summary, the prospects for the near-term first discoveries 
of Galactic TeV neutrino sources are very good. Importantly, 
this conclusion is empirically based on the measured spectra 
of bright Galactic TeV gamma-ray sources. For some of these, 
there are very compelling independent arguments that the ob- 
served gamma rays arise from neutral pion decays, meaning 
that they must be accompanied by neutrinos. It is essential 
to test this directly for these sources, as well as for others 
where the possibility of neutrino emission is uncertain. We 
emphasize the importance of using the measured muon en- 
ergy spectra to discriminate against the quickly falling atmo- 
spheric neutrino backgrounds. For example, a single event 
near 10 TeV from a source direction is almost certainly sig- 
nal, while an event near 1 TeV has a much higher probability 
of being background. This fact alone could be enough to help 
establish discovery. Due to the amplifying factors of neutrino 
cross section and muon range, neutrino detectors have better 
reach to the highest energies in the source spectra, as com- 
pared to gamma-ray telescopes, which can make precise mea- 
surements at lower energies. This complementarity can be 
exploited to help solve the long-standing puzzle of the origin 
of the Galactic cosmic rays. 

For all of the sources discussed here, the rates in km 3 
neutrino telescopes are relatively small, at most ~ 1 — 
10 events/year, though we have argued that even these small 
rates could have a powerful impact. While it remains possible 
that various uncertainties and detector limitations may make 
these observations even more challenging, they might also be 
better than described here. If the source spectra extrapolate 
to higher energies than we have assumed, or especially if the 
emitted TeV gamma-ray spectra have been diminished by ab- 
sorption in the sources, which is quite possible, then the neu- 
trino detection event rates could be significantly larger than 
shown here. Given the potentially unique power of neutrino 
astronomy, we can only hope that Nature has been so kind. 
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